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Using samples of 102 million T(IS') and 158 million T(2S) events collected with the Belle de- 
tector, we study exclusive hadronic decays of these two bottomonium resonances to the three- 
body final states <f>K + K~ , uj-k + -k~ and K* c '(892)K~ ty + ', and to the two-body Vector- Tensor 
states (0/2(1525), w/ 2 (1270), pa 2 (1320) and K*° (892) (1430)) and Axial-vector-Pseudoscalar 
(#1 (1270)+ if", ^1(1400) + ^" and &i(1235) + tt-) states. Signals are observed for the first time 
in the T(15) -> 4>K + K~ , u^ix" , K*°K~ir + , K*°Kf and T(25) ->■ 4>K + K~ , K*°K~ir + decay 
modes. Branching fractions are determined for all the processes, while 90% confidence level upper 
limits are established on the branching fractions for the modes with a statistical significance less 
than 3<r. The ratios of the branching fractions of T(2S) and Y(IS') decays into the same final state 
are used to test a perturbative QCD prediction for OZI suppressed bottomonium decays. 

PACS numbers: 13.25.Gv, 14.40.Pq, 12.38. Qk 



Although around 80% of the T(15) decays and 60% of 
the T(25) decays are expected to result in hadronic final 
states via annihilation into gluons [l], Q , no single exclu- 
sive mode has been reported |3J. This situation is quite 
different from the charmonium sector, where numerous 
channels have been measured and used to test a variety 
of theoretical models. The OZI (Okubo-Zweig-Iizuka) [4] 
suppressed decays of the J ftp and ip(2S) to hadrons pro- 
ceed via the annihilation of the quark-antiquark pair into 
three gluons or a photon. For both cases, perturbative 
quantum chromodynamics (pQCD) provides a relation 
for the ratios of branching fractions (B) for J/ip and 
ip(2S) decays [5J 
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which is referred to as the "12% rule" and is expected 
to apply with reasonable accuracy to both inclusive and 
exclusive decays. However, it is found to be severely vi- 
olated for pir and other Vector-Pseudoscalar (VP) and 



Vector- Tensor (VT) final states. None of the many ex- 
isting theoretical explanations that have been proposed 
have been able to accommodate all the measurements 
reported to date @. 

A similar rule can be derived for OZI-suppressed bot- 
tomonium decays, in which case we expect 
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0.77 ±0.07. 



T(1S)- 



(2) 

This rule should hold better than the 12% rule for char- 
monium decay, since the bottomonium states have higher 
mass, pQCD and the potential models should be more 
applicable, as has been demonstrated in calculations of 
the bb meson spectrum. 

In this Letter, we report studies of exclusive 
hadronic decays of the T(15) and T(2S) resonances 
to the three-body final states (f>K + K~, ujtt + tt~ , and 
K*°(892)K-tt + [71 and two-body VT [^(1525) -> 
0K+K-, cj/ 2 (_1270) ->• W7T+7T-, pa 2 (1320) -> p°p + n-, 
and if* (892)iY 2 *°(1430) -> AT*°(892)A'-7r+] and Axial- 
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vector-Pseudoscalar (AP) [^(1270)+^- p K + K~, 
Ki{U00) + K~ -> K*°(892)ir+K-, and &i(1235)+tt- 

final states. This analysis is based on a 
5.7 ft.- 1 T (IS) data sample (102 million T(15) events), a 
24.7 ftr 1 T(2S) data sample (158 million T(25) events), 
and a 89.4 fb^ 1 continuum data sample collected at 
y/s = 10.52 GeV. Here, y/s is the center-of-mass (CM.) 
energy of the colliding e + e~ system. The data are col- 
lected with the Belle detector [8( operating at the KEKB 
asymmetric-energy e + e~ collider [9j. The evtgen [lOj 
generator is used to simulate Monte Carlo (MC) events. 
For two-body decays, the angular distributions are gen- 
erated using the formulae in Ref. ll|. Inclusive T(15) 
and T(2S) MC events, produced using PYTHIA [ijj with 
the same luminosity as real data, are used to check for 
possible peaking backgrounds. 

We require four reconstructed charged tracks with zero 
net charge. For these tracks, the impact parameters 
perpendicular to and along the beam direction with re- 
spect to the interaction point are required to be less than 
0.5 cm and 4 cm, respectively, and the transverse momen- 
tum in the laboratory frame is restricted to be higher 
than 0.1 GeV/c. For each charged track, we combine in- 
formation from different detector subsystems to form a 
likelihood Ci for each particle species [131 ] . 

A track with TZk = Cr+c ^ ^-6 is identified as a 
kaon, while a track with TZk < 0.4 is treated as a pion. 
With this selection, the kaon (pion) identification effi- 
ciency is about 89% (92%), while 6% (9%) of kaons (pi- 
ous) are misidentified as pions (kaons). A similar likeli- 
hood ratio Hp is formed for muon identification [l4| . Ex- 
cept for the (f>K + K~ final state for which at least three 
charged tracks are required to be identified as kaons, all 
other charged tracks are required to be positively identi- 
fied as pions or kaons. A small background with muons is 
removed by requiring 71^ < 0.95 for the pion candidates. 

A neutral cluster in the electromagnetic calorimeter 
is reconstructed as a photon, if it does not match the 
extrapolated position of any charged track and its energy 
is greater than 40 MeV. A tt° candidate is reconstructed 
from a pair of photons. We perform a mass-constrained 
fit to the selected ir° candidate and require x 2 < 6. 

We impose an energy conservation requirement on 
Xt = Y^hEh/ \fs, where Eh is the energy of the final 
state particle h in the e + e~ CM. frame. The ratio 
X T should lie in the range 0.985 < X T < 1.015 for 
channels with a ir° in the final state, and in the range 
0.99 < Xt < 1.01 for other channels. Figure [T] shows 
the X T distributions from T(15) and T(25) decays to 
0K+K-, ujtt+tt- and if*°(892)iC7r + , together with ex- 
pected backgrounds from continuum processes and Y(1S) 
and Y(2S) decays. Obvious signal candidates at Xt ~ 1 
can be seen. 

For three-body decay modes with a (f> or w, <j) (w) 
candidates are selected with K + K~ (7r + 7r _ 7r°) masses 
closest to the nominal <f> (oj) mass [![. Figure [5] shows 



30 
20 
> 10 





LU 



T 



Y(1S)^*K*K" 



• r 


| y(2S)^*k*k: 






J 








4oJ 



0.96 0.98 



1.02 1.04 



0.96 0.98 



1 



1 .02 1 .04 



£ 10 















I|y(2S) -> K'°K"tt*: 






if : 















0.96 0.98 



FIG. 1: Scaled total energy, Xt, distributions from T(1S) 
and T(25) decays to <j>K + K~, wtt+tt" and K*° (892)K~ ir+ . 
The dots with error bars are from resonance data; the dark 
shaded histograms are from normalized continuum contribu- 
tions; the light shaded histograms are from inclusive T(1S) 
and T(2S) MC events with signals removed. The arrows show 
the required signal region. 



the K + K , ir + ir 7r°, and K + ir invariant mass distri- 
butions for T(1S) and T(2S) to 4>K+K-, ujtt+it- and 
K*°(892)K ~7r + candidates that survive the selection cri- 
teria described above. Clear cj), u, and if* (892) signals 
are evident. 

After the application of all of the selection require- 
ments, no peaking backgrounds from the T(1S) and 
T(2S) inclusive MC samples are found in the vector me- 
son mass regions. Potential backgrounds due to parti- 
cle misidcntification, from (/>7t + 7t - for example, are esti- 
mated by selecting these events in the data and normal- 
izing them using measured misidentification probabili- 
ties. Potential backgrounds from events with additional 
7r°'s are checked by examining the recoil mass distribu- 
tion from the measured final state. For W7r + 7r~ , potential 
background events from iorj with if — > 77r + 7r~ are explic- 
itly reconstructed from data and estimated using Nei/e 2 , 
where N is the number of ojr]' events in data and e± and 
62 are the efficiencies after the W7r + 7r~ and uirf event se- 
lections, respectively. All of the above backgrounds are 
found to be negligibly small. For um + ir~, the fraction 
of events with multiple combinations is at the 3.5% level 
due to multiple ir° candidates; this is consistent with the 
MC simulation and is taken into account in the efficiency 
determination. 
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The continuum background contribution is determined 
from the data at y/s = 10.52 GeV and is extrapo- 
lated down to the T(IS') and T(25) resonances. For 
the extrapolation, the scale factor, / sca i e j is given by 
Ct — ctt — cr_ where Or_ an d are the ratios 

t con (T con ^con L con (T con ' €con 

of luminosity, cross sections, and efficiencies at the bot- 
tomonium masses and continuum energy points. The s 
dependence of the cross section is assumed to be 1/s [HI 
and the corresponding scale factor is 0.079 for the T(15) 
and 0.30 for the T(25). 

An unbinned simultaneous likelihood fit to the mass 
distributions is applied to extract the signal and back- 
ground yields in the T(IS') and continuum data samples 
and in the Y(2S) and continuum data samples. The sig- 
nal shapes are obtained from MC simulations. In this fit, 
a second-order Chebyshev polynomial background shape 
is used for the T(1S)/T(2S) decay backgrounds in addi- 
tion to the normalized continuum contribution. The fit 
ranges and results for the K + K~, Tr + ir~ir°, and K + ir~ 
mass spectra are shown in Fig. [2] and Table Q] 

We determine a Bayesian 90% confidence level (C.L.) 
upper limit on N sl& by finding the value such that 

/°° cdN si s = 0-90, where iV s i g is the number of signal 
events and C is the value of the likelihood as a function 
of iV S ig. The statistical significance of the signal is esti- 
mated from the difference of the logarithmic likelihoods, 
—2 ln(£o/£ ma x)j taking into account the difference in the 
number of degrees of freedom in the fits, where Cq and 
£ max are the likelihoods of the fits with and without sig- 
nal, respectively. 

After requiring \M K+K - — m,p\ < 8 MeV/c 2 , 
\M„+ v -„o - m u \ < 30 MeV/c 2 , and \M K+n - - 
w_R"*°(892)| < 100 MeV/c 2 , which contain around 95% 
of the signal according to MC simulations, the Dalitz 



plots for the (j)K + K~ 



and K*°(892)K-ir + final 



states are shown in Fig. [3] where m,p, m^, and m^*o( 8g2 ) 
are the nominal <f>, w, and K* a (892) masses [l[. Interest- 
ingly, the events accumulate near the phase space bound- 
ary, reflecting the quasi-two-body nature of these decays. 

To categorize the quasi-two-body decays into VT or 
AP final states, we further require the angle between V 
and T (A and P) in the e + e~ CM. frame to be greater 
than 179 degrees for the channels with a 7r°, or 179.5 
degrees for the other channels. The combination with 
the minimum value of 8 m i n = {Mi — my) 2 + (M% — mx) 2 
is selected as the V and T candidate, where Mi and M 2 
are the invariant masses of the V and T decay final-state 
particles, respectively. The same technique is used to 
select the best AP candidate. This method introduces 
negligible bias in the meson pair selection according to 
MC simulation. 

For the selected events, Fig.|4]shows the invariant mass 
distributions for the vector and tensor meson candidates, 
and Fig. [3] shows the invariant mass distributions for the 
vector (from axial- vector decay) and axial- vector meson 
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FIG. 2: The K + K~ (top row), tt+tt^tt (middle row), and 
K + tt~ (bottom row) invariant-mass distributions for the final 
candidate events from Y(IS') (left column) and T(2S) (right 
column) three-body decays. Solid points with error bars are 
data, open histograms show the best fits, dashed curves are 
the total background estimates, and shaded histograms are 
the normalized continuum background contributions. 
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FIG. 3: Dalitz plots of 4>K + K~ (top row), uj-k + tt~ (mid- 
dle row) and K*° (892)K~ ir + (bottom row) three-body final 
states. Here, the left column is for Y(IS') decays, the middle 
column is for T(2S) decays, and the right column is for the 
continuum data. 
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TABLE I: Results for the T(1S) and T(2S) decays, where N B1B is the number of fitted signal events, is the upper limit 

on the number of signal events, e is the efficiency (%), E is the statistical significance (a), B is the branching fraction, B VP is 
the upper limit on the branching fraction, Qr is the ratio of the T(2S) and T(15) branching fractions, and Q^ p is the upper 
limit on the value of Qr- Branching fractions are in units of 10 -6 and upper limits are given at the 90% C.L. The first error 
in B and Qr is statistical, and the second systematic. 



Channel 


T(1S) 

JV sig e £ B B up 


T(2S) 

N Bie nV? t £ B B l,F 


Qr Q? P 


4>K+K~ 
K'°K-k+ 


56.3 ± 8.7 47.9 8.6 2.36 ± 0.37 ± 0.29 
63.6 ± 9.5 15.7 8.5 4.46 ± 0.67 ± 0.72 
173 ± 20 28.7 11 4.42 ± 0.50 ± 0.58 


58 ± 12 47.8 6.5 1.58 ± 0.33 ± 0.18 

29 ± 12 51 15.9 2.5 1.32 ± 0.54 ± 0.45 2.58 

135 ± 23 27.5 6.4 2.32 ± 0.40 ± 0.54 


0.67 ± 0.18 ± 0.11 

0.30 ± 0.13 ± 0.11 0.55 

0.52 ± 0.11 ± 0.14 


<f>fL 

pa.2 
K*°K* 


6.9 ± 3.9 15 48.8 2.1 0.64 ±0.37 ±0.14 1.63 
5.2 ± 4.0 13 17.7 1.5 0.57 ± 0.44 ± 0.13 1.79 
29 ± 11 49 17.4 2.7 1.15 ± 0.47 ± 0.18 2.24 
42.2 ± 9.5 30.8 5.4 3.02 ± 0.68 ± 0.34 


8.3 ±6.0 18 49.0 1.6 0.50 ± 0.36 ± 0.19 1.33 
-0.4 ± 3. 3 6.1 17.5 -0.03 ± 0.24 ± 0.01 0.57 
10 ± 11 30 17.3 0.9 0.27 ± 0.28 ± 0.14 0.88 
32 ± 11 29.6 3.3 1.53 ± 0.52 ± 0.19 


0.77 ± 0.70 ± 0.33 2.54 
-0.06 ± 0.42 ± 0.02 1.22 
0.23 ± 0.26 ± 0.12 0.82 
0.50 ± 0.21 ± 0.07 


ifi(1270) + if~ 
Ki(Um) + K- 
hi(1235) + 7r~ 


3.7 ± 4.9 13 23.6 0.8 0.54 ± 0.72 ± 0.21 2.41 
23.8 ± 8.2 27.3 3.3 1.02 ± 0.35 ± 0.22 
14.4 ± 6.9 28 16.7 2.4 0.47 ± 0.22 ± 0.13 1.25 


11.0 ± 4.4 26 23.5 1.2 1.06 ± 0.42 ± 0.32 3.22 
9.2 ± 8.2 24 26.9 0.5 0.26 ± 0.23 ± 0.09 0.83 
1.2 ± 3.5 13 17.0 0.2 0.02 ± 0.07 ± 0.01 0.40 


1.96 ± 2.71 ± 0.84 4.73 
0.26 ± 0.25 ± 0.10 0.77 
0.05 ± 0.16 ± 0.03 0.35 



candidates for T(IS') and T(2S) two-body decays. We 
extend the unbinned simultaneous maximum likelihood 
fit described above for three-body decays into a two- 
dimensional (2D) fit. We assume the mass distributions 
of the V and T particles to be uncorrelated; thus, the 
mass distributions in the 2D space can be represented 
by the product of two one-dimensional (ID) probability 
density functions (pdf). The 2D fitting function is pa- 
rameterized as 

f(M lt M 2 ) = iV si s Sl (Af 1 )s 2 (M 2 ) + N h ^s l {M l )b 2 {M 2 ) 
+N^b 1 (M 1 )s 2 {M 2 ) + iV b b h s 6i(Afi)6 2 (M 2 ), 

where s\(Mi) and bi(M\) are the ID signal and back- 
ground pdfs for V, respectively, and s 2 (M 2 ) and b 2 (M 2 ) 
are the corresponding pdfs for T . Here, the free parame- 
ters are the signal yield N S1S and the background yields 
N sb> N bs> and N bb- Similar 2D pdfs are used to fit the 
vector (in axial-vector decays) and axial-vector meson 
candidates for the AP modes. In these fits, we assume 
there is no interference between the signal and other com- 
ponents due to the limited statistics. The ID projections 
from the 2D fits are shown in Figs. [4] and [5] with the con- 
tribution from each component indicated. In the fits to 
the #i(1270)+if- and Ki(1400)+X" modes, the cross- 
feed background components from ifi(1400) + K~ — » 
K*°K-tt+ and ^(1270)+ K~ -> p°K+K~ are also in- 
cluded and shown as dot-dashed lines. The fit results are 
shown in Table HI 

There are several sources of systematic errors for the 
branching fraction measurements. The uncertainty in the 
tracking efficiency for tracks with angles and momenta 
characteristic of signal events is about 0.35% per track 
and is additive. The uncertainty due to particle identifi- 
cation efficiency is 1% with an efficiency correction factor 
of 0.97 for each pion and 0.8% with an efficiency correc- 
tion factor 0.97 for each kaon, respectively. The uncer- 
tainty in selecting 7r° candidates is estimated by compar- 



ing control samples of rj — > 7r 7r 7r° and r\ — > 7r + 7r~7r° de- 
cays in data and amounts to 3.7%. Errors on the branch- 
ing fractions of the intermediate states are taken from 
the PDG listings According to MC simulation, the 
trigger efficiency is greater than 99%, so that the corre- 
sponding uncertainty can be neglected. We estimate the 
systematic errors associated with the fitting procedure 
by changing the order of the background polynomial and 
the range of the fit; the differences in the fitted results, 
which are 1.3%-29% depending on the final state parti- 
cles that are taken as systematic errors. We estimate the 
systematic errors associated with the resonance param- 
eters by changing the values of the masses and widths 
of the resonances by ±Icr; the differences of 0.6%-7.3% 
in the fitted results are taken as systematic errors. For 
the central values of the branching fractions, the average 
difference between alternative CM. energy dependences 
of the cross section is included as a systematic error due 
to the uncertainty of the continuum contribution, which 
is in the range of 4.2% to 22%. The uncertainty due to 
limited MC statistics is at most 0.5%. Finally, the uncer- 
tainties on the total numbers of T(f S) and T(2S) events 
are 2.2% and 2.3%, respectively. Assuming that all of 
these systematic error sources are independent, the total 
systematic error is 11%-31% depending on the final state. 

Table Q] shows the results for the branching fractions 
including the upper limits at 90% C.L. for the channels 
with a statistical significance less than 3er. The corre- 
sponding ratio of the branching fractions of T(2S) and 
T(IS') decay is calculated; in some cases, the systematic 
errors cancel. In order to set conservative upper limits 
on these branching fractions, the efficiencies are lowered 
by a factor of I — <7 sys in the calculation, where er sys is 
the total systematic error. All the results on the branch- 
ing fractions, including upper limits, are below CLEO's 
preliminary results 

In summary, we have measured T(IS') and T(2S) 
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FIG. 4: The mass projections for the vector and tensor meson candidates from 2D fits to the events from T(IS') and T(2S) 
two-body decays (VT modes). The open histograms show the results of the 2D simultaneous fits, the dotted curves show the 
total background estimates, and the shaded histograms are the normalized continuum contributions. 
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FIG. 5: The mass projections for the vector (from axial-vector decay) and axial-vector meson candidates from 2D fits to the 
events from T(IS') and T(2S) two-body decays (AP modes). The open histograms show the results of the 2D simultaneous 
fits, the dotted curves show the total background estimates, the dot-dashed curves are the cross-feed backgrounds described in 
the text, and the shaded histograms are the normalized continuum contributions. 
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hadronic exclusive decays to three-body final states and 
two-body processes. Signals are observed for the first 
time in the T(15) -> 4>K+R-, cjtt+tt-, K*°K-tt+, 
K*°Kf and T(25) -> (j)K + K~, K*°K-tt+ decay 
modes. Besides iC* ^ , no other two-body processes are 
observed in all investigated final states. We find that for 
the processes cj>K+K-, K*°K-ir+ 1 and K*°Kf(U30), 
the Qt ratios are consistent with the expected value, 
while for W7r + 7r~, the measured Qx ratio is 2.6a below 
the pQCD expectation. The results for the other modes 
are inconclusive due to low statistical significance. These 
results may supply useful guidance for interpreting vio- 
lations of the 12% rule for OZI suppressed decays in the 
charmonium sector. 
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